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Abstract 

The experimental charged particles multiplicity distributions in e^e^ anni- 
hilation to hadrons are compared with the distributions obtained by PYTHIA. 
The ratio x^/degrees of freedom is calculated for 6 energies at ^/s 14 - 206.2 
GeV. The necessity of more subtle tuning of PYTHIA at the energy of peak 
is discussed. 

1 Introduction 

In the last decades both experimentalists and theorists have applied significant forces 
to understand character of hadron interactions at high energies. As result many as- 
pects of hadron processes at high energies became clear not only qualitatively but also 
quantitatively. The last corresponds to quantum chromodynamics which was able to 
describe quantitatively many the simplest processes occurring at short distances. How- 
ever, it is impossible to count on QCD calculations of the multiple hadron production, 
these processes have enormous number of degrees of freedom. Thus one has to use 
different model approaches. A quantitative description of high energy processes by 
Monte Carlo generators can be considered as one of these approaches. Monte Carlo 
generators are perspective method of theoretic study at all energies of the existing ac- 
celerators including LHC, since any other approaches, theoretical as the BFKL theory 
or phenomenological as the Regge theory, are only valid to the values of 1/ In ^/s order 
of magnitude, where s is the square of total center-of-mass energy. Since even for LHC 
In -^i ~ 19, then it is necessarily to take into account the next to 1/ In ^/s corrections 
what is an inexecutable task at the most cases. Therefore it is essentially to construct 
such MC generator which will describe well processes both in hard and soft regions. 

Amongst generators describing well hard processes there is known PYTHIA [I]. Still 
soft processes, i.e. processes with small transverse momenta, are included in PYTHIA 
only on basis of the simplest phenomenological models, the values of cross sections are 
not generated. The multiparticle processes are described with the simplest pomeron 
model in which only one string is stretched between quark and diquark. Pomeron 
cuts are not taken into account. The weakness of PYTHIA is also in absence of the 
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AGK sum rules [2]. An oscillation structure in the multiplicity distribution function 
Pn = <^n/o'tot, where cr„ is the topological cross section of n hadrons production, dtot is 
the total cross section of hadron-hadron scattering, that can be observed at the LHC 
energies, follows from the AGK rules. Since PYTHIA is the most advanced programm 
then the future MC generator describing both hard and soft regions, in our opinion, 
should be based on PYTHIA. 

In the present work we will study how multiplicity distributions obtained with 
PYTHIA agree with experimental data. 

2 The charged particles multiphcity distribution in 
e+e~ annihilation to hadrons and PYTHIA pre- 
dictions 

We consider the process of e~^e~ annihailation at 6 values of total center-of-mass energy 
from 14 to 206.2 GeV. Wide coverage of the experimental data and sufficient statistical 
provision of them were criteria for the choice. 

Electron-positron annihilation to hadrons was simulated by PYTHIA at ^/s 14, 29, 
34.8, 91.2, 188.6, 206.2 GeV. The decay process of 7* (Z") was considered without the 
initial state radiation, values of the other parameters were set by default. There were 
generated one million events for every energy, and the charged particles multiplicity 
distributions were analyzed. 
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2.1 v^=14GeV 



Experimentally there were obtained 2 704 events (TASSO) [3]. The experimental values 
of Pn and the corresponding probabilities obtained by PYTHIA are shown in Fig.l. 
Both statistical and systematical errors are included. The is presented in Table 1. 
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Table 1. y/i = 14 GeV 



n 


number of events, experiment 


number of events, PYTHIA 


in bin 


0-2 


15.64 


12.20 


0.33 


4 


145.98 


116.08 


3.43 


6 


445.47 


414.55 


1.20 


8 


720.66 


720.78 


0.00 


10 


677.63 


707.34 


0.71 


12 


405.05 


443.66 


2.06 


14 


181.75 


194.74 


0.55 


16 


70.79 


66.74 


0.15 


18 


23.53 


20.51 


0.28 


20 - 30 


17.51 


7.40 


2.64 


total x^/degrees of freedom = 11/9 
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2.2 = 29 GeV 



Experimentally there were obtained 29 649 events (HRS) [4J. The experimental values 
of Pn and the corresponding probabilities obtained by PYTHIA are shown in Fig.2. 
Both statistical and systematical errors are included. The is presented in Table 2. 
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Table 2. = 29 GeV 



n 


number of events, experiment 


number of events, PYTHIA 


in bin 


0-2 


14.82 


12.87 


0.05 


4 


186.79 


169.59 


0.08 


6 


1055.50 


952.57 


1.98 


8 


2795.90 


2811.24 


0.01 


10 


5200.43 


5122.86 


0.10 


12 


6854.85 


6313.39 


2.91 


14 


5760.80 


5724.38 


0.02 


16 


3949.25 


4057.80 


0.32 


18 


2188.10 


2397.37 


2.73 


20 


1064.40 


1212.53 


2.82 


22 


367.65 


545.37 


11.70 


24 


148.25 


213.92 


4.70 


26 


50.40 


79.13 


2.76 


28 - 36 


11.86 


35.99 


7.29 


total x^/degrees o 


■ freedom = 37/13 
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2.3 = 34.8 GeV 



Experimentally there were obtained 52 832 events (TASSO) [3J. The experimental 
values of P„ and the corresponding probabilities obtained by PYTHIA are shown in 
Fig. 3. Both statistical and systematical errors are included. The is presented in 
Table 3. 
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Table 3. y/i = 34.8 GeV 



n 


number of events, experiment 


number of events, PYTHIA 


in bin 


0-2 


23.62 


22.93 


0.00 


4 


302.89 


302.20 


0.00 


6 


1673.45 


1697.39 


0.08 


8 


4427.16 


5009.39 


19.14 


10 


8090.16 


9128.50 


33.12 


12 


10456.88 


11249.92 


14.78 


14 


10167.52 


10200.36 


0.03 


16 


7625.51 


7230.66 


5.23 


18 


4850.98 


4271.90 


18.46 


20 


2674.51 


2160.62 


27.65 


22 


1434.97 


971.79 


43.62 


24 


652.37 


381.18 


33.72 


26 


273.30 


141.01 


19.60 


28 


104.45 


47.23 


9.52 


30 - 36 


74.28 


16.91 


4.48 


total /degrees of freedom = 229/14 
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2.4 v^ = 91.2GeV 



Experimentally there were obtained 248100 events (L3) [5J. The experimental values 
of Pn and the corresponding probabilities obtained by PYTHIA are shown in Fig. 4a. 
Both statistical and systematical errors are included. The is presented in Table 4a. 

Unlike the previous ones, the ratio of x^/degrees of freedom = 1155/25 is very large 
for the results of L3 col., although visually the experimental and generated points are 
in good agreement (especially if one takes the curve fitting the generated points but 
not the points themselves, as it was done in some articles, for example [S], [B]). This 
is bound to very large statistics in the L3 col. and its small errors. The multiplicity 
distributions for the same energy obtained in DELPHI [7J (25 364 events) and OPAL [8] 
(82 941 events) experiments are given below (Fig.4b, Table 4b and Fig.4c, Table 4c 
accordingly). In these experiments statistics is lower and the values of x^/degrees of 
freedom are correspondingly smaller. The result of multiplicity distributions from all 
collaborations combined together is given in Fig.4d. 
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Table 4a. ^ = 91.2 GeV, L3 col. 



n 


number of events, experiment 


number of events, PYTHIA 


in bin 


2-4 


70.96 


39.94 


21.15 


6 


509.60 


343.37 


54.22 


8 


2314.28 


1602.97 


218.37 


10 


6852.77 


5177.60 


203.97 


12 


14495.49 


11806.33 


151.89 


14 


23280.71 


20305.99 


90.10 


16 


30198.48 


27721.70 


48.02 


18 


33190.57 


32102.16 


8.85 


20 


32138.63 


32137.31 


0.96 


22 


28183.66 


29575.26 


24.48 


24 


22970.59 


24935.54 


104.10 


26 


17573.42 


19586.50 


88.30 


28 


12708.92 


14481.85 


56.91 


30 


8805.07 


10245.29 


28.70 


32 


5840.52 


6955.73 


21.36 


34 


3703.64 


4480.69 


13.01 


36 


2280.04 


2734.06 


7.44 


38 


1338.50 


1619.84 


5.33 


40 


773.82 


920.45 


3.18 


42 


430.45 


514.31 


1.97 


44 


224.03 


258.77 


0.76 


46 


112.14 


134.97 


0.81 


48 


56.57 


68.23 


0.51 


50 


25.31 


31.51 


0.35 


52 


12.65 


12.90 


0.00 


51 62 


8.93 


6.70 


0.07 


total x^/degrees of freedom = 1155/25 
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Table 4b. = 91.2 GeV, DELPHI col. 



n 


number of events, experiment 


number of events, PYTHIA 


in bin 


2-6 


45.91 


39.19 


0.23 


8 


170.95 


163.88 


0.14 


10 


578.30 


529.32 


1.10 


12 


1230.15 


1207.00 


0.09 


14 


2084.92 


2075.94 


0.01 


16 


2815.40 


2834.07 


0.01 


18 


3271.96 


3281.90 


0.00 


20 


3322.68 


3316.17 


0.00 


22 


2967.59 


3023.57 


0.12 


24 


2483.14 


2549.23 


0.23 


26 


1909.91 


2002.39 


0.70 


28 


1460.97 


1480.52 


0.05 


30 


1050.07 


1047.41 


0.00 


32 


743.17 


711.11 


0.40 


34 


476.84 


458.07 


0.28 


36 


309.44 


279.51 


1.30 


38 


191.50 


165.60 


1.83 


40 


121.24 


94.10 


1.00 


42 


63.66 


52.58 


0.43 


44 


36.27 


26.45 


0.92 


46 


20.80 


13.80 


1.16 


48 


5.07 


6.98 


0.39 


50-62 


4.31 


5.22 


0.02 


total x^/ degrees o 


f freedom = 10/22 
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Table 4c. ^ = 91.2 GeV, OPAL col. 



n 


number of events, experiment 


number of events, PYTHIA 


in bin 


2-4 


14.10 


13.35 


0.00 


6 


132.71 


114.79 


0.04 


8 


564.00 


535.88 


0.03 


10 


1725.17 


1730.90 


0.00 


12 


3889.93 


3946.91 


0.06 


14 


6635.28 


6788.39 


0.48 


16 


8949.33 


9267.50 


0.71 


18 


10458.86 


10731.90 


0.99 


20 


10657.92 


10813.96 


0.29 


22 


9811.92 


9887.15 


0.09 


24 


8285.81 


8336.07 


0.02 


26 


6510.87 


6547.86 


0.04 


28 


4934.99 


4841.35 


0.24 


30 


3607.93 


3425.05 


1.06 


32 


2463.35 


2325.33 


1.41 


34 


1675.41 


1497.91 


2.75 


36 


1069.94 


914.01 


2.08 


38 


671.82 


541.52 


3.58 


40 


389.82 


307.71 


2.16 


42 


215.65 


171.94 


1.01 


44 


141.00 


86.51 


1.43 


46 


73.82 


45.12 


0.68 


48 


34.84 


22.81 


0.38 


50 


20.74 


10.53 


0.47 


52 - 62 


12.44 


6.55 


0.20 


total \;^/(legrees o 


f frecxlom = 20/21 
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2.5 = 188.6 GeV 



Experimentally there were obtained 4479 events (L3) [H]. The experimental values of 
Pn and the corresponding probabilities obtained by PYTHIA are shown in Fig.5. Both 
statistical and systematical errors are included. The is presented in Table 5. 
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Table 5. ^ 188.6 GeV 



n 


number of events, experiment 


number of events, PYTHIA 


in bin 


2-10 


17.92 


18.16 


0.00 


12 


58.23 


44.23 


2.32 


14 


125.41 


101.23 


2.23 


16 


214.99 


186.12 


1.87 


18 


322.49 


279.53 


3.41 


20 


389.67 


364.89 


1.09 


22 


434.46 


425.13 


0.11 


24 


452.38 


447.01 


0.04 


26 


131.16 


111.80 


0.05 


28 


394.15 


413.00 


0.47 


30 


353.84 


366.51 


0.21 


32 


309.05 


317.03 


0.09 


34 


255.30 


260.58 


0.03 


36 


201.56 


210.87 


0.11 


38 


152.29 


165.96 


0.25 


40 


116.45 


125.96 


0.13 


42 


80.62 


93.32 


0.37 


44 


58.23 


68.70 


0.27 


46 


40.31 


48.67 


0.18 


48 


22.40 


33.54 


2.32 


50 


17.92 


23.68 


0.09 


52 - 82 


8.96 


43.07 


8.12 


total x^/degrees o 


f freedom = 24/21 
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2.6 = 206.2 GeV 



Experimentally there were obtained 4146 events (L3) [H]. The experimental values of 
Pn and the corresponding probabilities obtained by PYTHIA are shown in Fig.6. Both 
statistical and systematical errors are included. The is presented in Table 6. 



0.1- 



0.08 



0.06 ^ 



0.04 



0.02 



.it 



— i-DGa-n-i 

J_l I I I I I I L_L 



206.2 GeV, x2/dof=1 7/21 
• Experiment 
□ Pythia 



J I I I I I I LJ I Ll I I I I I I LJ L 



J_J LJ I I LJ I Ll I I I 



10 20 30 40 50 



60 70 80 90 

n 



Fig.6. 



15 



Table 6. ^ 206.2 GeV 



n 


number of events, experiment 


number of events, PYTHIA 


in bin 


2-10 


12.44 


13.43 


0.02 


12 


41.46 


32.33 


1.25 


14 


91.21 


77.09 


1.22 


16 


165.84 


146.40 


1.02 


18 


257.05 


229.01 


1.74 


20 


344.12 


306.95 


2.89 


22 


385.58 


366.97 


0.53 


24 


402.16 


395.36 


0.07 


26 


389.72 


100.01 


0.11 


28 


364.85 


381.23 


0.31 


30 


319.24 


347.40 


1.00 


32 


281.93 


303.93 


0.60 


34 


240.47 


259.44 


0.51 


36 


194.86 


214.66 


0.55 


38 


157.55 


172.92 


0.35 


40 


128.53 


135.30 


0.07 


42 


99.50 


103.13 


0.02 


44 


74.63 


78.00 


0.03 


46 


58.04 


56.78 


0.00 


48 


49.75 


40.26 


0.42 


50 


33.17 


28.50 


0.05 


52-88 


24.88 


56.91 


4.48 


total x^/degrees o 


f freedom = 24/21 
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3 Conclusion 



As it was shown, PYTHIA describes well the experimental data on the multiplicity 
distribution in e^e~ annihilation to hadrons. So it can be used as a reliable method of 
theoretical studies. 

As it was noticed in PYTHIA manual [1], the default values of all its parameters 
were tuned to describe the events at the energy of peak. The ratio of x^/degrees 
of freedom is equal to 10/22 for DELPHI col. with 25 364 events, 20/24 for OPAL 
col. with 82 941 events and it is surprisingly large, 1155/25, for L3 col. with 248100 
events. The ratio is approximately equal to 1 for other discussed energies, although 
for 34.8 it is 229/14 (52 832 events, TASSO col.). This behavior of is quite obvious. 
Thus, more subtle tuning is necessary for the values of PYTHIA parameters, so that 
PYTHIA will be able to predict quantity results. 

Author is grateful for V. A. Abramovsky for useful discussions. 
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